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We obtain analytical solutions for positive-to-negative transition materials with a 
smooth transition profile "tanh" . The fields are expressed in terms of hypergeomet- 
ric functions. The final expressions for reflection and transmission coefficients are 
obtained in closed forms. The total absorption is nonzero even in the lossless limit in 
accordance with previous studies. Properties of the total absorption as function of 
controlling factors are also well studied. As an application of the analytical results, 
we analyze the effects of smooth transition on perfect imaging and, interestingly, find 
multiple images of the object. 



I. INTRODUCTION 



Positive-to-negative transition materials are layered optical materials with gradient opti- 
cal indexes e and/or fi continuously changing from positive values to negative ones. In Refill 
and 2, it is obtained that near the transition point for e, the fields for oblique incidence 
present large enhancement with enhanced absorption. Actually the absorption is nonzero 
even in the lossless limit. In Ref. [j, similar effects are obtained when both e and p linearly 
pass the same transition point, and the fields are expressed in terms of confluent hypergeo- 
metric functions near that point. Also in Ref. [j, smooth transition material e = /i = tanhx 
(tanh-material) is only studied numerically. In Ref. |4|, the properties for the "tanh" profile 
under normal incidence are studied, but no field enhancement is present for normal inci- 
dence. In this study, we obtain analytical solutions for the "tanh" material under oblique 
incidence. The fields are expressed in terms of hypergeometric functions, and the reflection 
as well as the transmission coefficients are expressed in closed-form analytical formulas. The 
absorption in the lossless limit is nonzero, as expected. Therefore, this study provides an 
example of transition materials^ that can be treated fully analytically 

As an application of the obtained formulas, we also analyze the effects of smooth transition 
on perfect imaging. In the sharp transition limit, the material can have perfect imaging 
effect, i.e., all evanescent waves can be ideally recovered^. However, for smooth transition, 
it is not that case. Our analytical results can shed light on the effects of smooth transition. 
The main feature of the image plane is that it contains multiple images of an object. 



II. HELMHOLTZ EQUATION FOR STEADY STATES 

The permittivity and permeability we are interested in are e = p = tanh(px + i5). The 
imaginary part is positive for < 5 < tt/2, and therefore the material is passive under that 
condition. 

The Helmholtz equation for TE waves are j^], 
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-— - -Jt— + aV - b 2 )E = 0. (1) 
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where a = u/ p, b = usinO/p with u and 9 being the frequency and angle of the incident 
wave. 

The Eqfl] can be transformed into hypergeometric equation after changes of variables, 



y = ji 2 and E — (1 — y) x F(y) where A = — ^\Ja? — b 2 . The resulted equation is 

F = 0, (2) 
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dy 2 l-y dy 4y(l - y) 
which is a hypergeometric equation with indexes a = A + ia/2, (3 = A — ia/2 and 7 = 0. 

III. SOLUTIONS FOR x -> +oo, x -> -oo AND -oo < x < +oo. 

For x — )■ ±oo, we have /x — )■ ±1 and y — > 1. According to Appendix A, the field solution 
can be expressed as linear superposition of two basic solutions, 

E ± {x) = - n 2 ) x w 3 + C¥\l - fi 2 ) V, (3) 

where iy 3j4 , as well as Wi )2 below, are Kummer solutions for hypergeometric equations in 
Appendix A. 

For — oo < x < +oo, y is near the zero pole. The two basic solutions are, 

(1 - M 2 ) A / and (l-/i 2 ) A w 2 - (4) 



where W2 is one of the Kummer solutions near the zero pole [10| and / is the solution with 
logarithm since we are working with a degenerate hypergeometric equation with index 7 = 0. 

Similar to the expression of Neumann functions in terms of Bessel function, the logarithm 
solution / can be expressed as limit of linear superposition of w\ and W2 of nonzero 7 index, 

- a/3w 2 

which will be proved in the Appendix B. Specifically, according to the L 'Hospital rule, 

Sr(7)-^ir'(7) _ a . 
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r( 7 )2 

where w[ = ^, w' 2 = ^ = -w 2 \ny + y 1 "^ 2 F 1 (a - 7 + 1, P - 7 + 1; 2 - T, v) with l F 2 
being the hypergeometric function. 

However, if we use the principal values of f(y), the field is discontinuous^ along the 
real x axis at zero point x = 0, corresponding to /x = tanh(z'5) = it&nhS. Therefore, for 
< x < +00 and —00 < x < 0, we should use different superpositions of w 2 and / to make 
sure that the field is continuous. 
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The discontinuity of the principal value of / is f(x = 0+) — f(x = 0— ) = 2a(3niw2(x = 0) 
according to Eq.(jHJ). Then, if the field for < x < +oo is 

E + (x) = A(l - fi 2 ) x f + B(l - /i 2 )V, (7) 

the field in — oo < x < should be 

E-(x) = A{1 - fi 2 ) x f + (B + 2af3rnA)(l - fi 2 ) x w 2 , (8) 

to ensure the continuity across x — 0. 



IV. RELATIONS AMONG C [ l' 2) AND A, B 



The different field expressions (0 [7f8|) should be equal in equivalent regions, which imposes 
connection relations among the coefficients C± ' and A, B. 

If we set = 1, = in x — > +oo region, according to the Kummer connection 
formula between W3 and w\ t 2 in Appendix A, as well as Eq. ( I51P7I) . we have 



w 3 



r(i - 7 )r(« + p - 7 + 1) , r( 7 -i)r(a + /3- 7 + i) 



r(« - 7 + l)T(/3 - 7 + 1) 
[5(7) - A( 7 )a/3r(7)]u; 2 + A( 7 ) Wl . 



r(«)r(/3) 



(9) 



with B = lim^o B(j) : A = lim^o M^i)- After straightforward calculation, we have 

T(a + f3 + 1) 



A 
B 



apV{a)V{p) 
1 



(- + ^ + 7J5 + #*) + V'G0)-l)- 



(10) 



T(a)T(l3) x a (3 

where je = — ^(1) = 0.577215664.... and ^ is the derivative of the logarithm of T function 
Then according to Eq. ( 171IH1I91) . in x — > — oo 



E-(x) = (1 - /i 2 ) A ^ 3 + 2a/37riA(l - /i 2 ) A w 2 



2\A„ 



(11) 



By use of Kummer connection formula between w<i and w 3 4 in Appendix A, we have 



Ejx) 



_ Aa/?(27ri)r(-a - /9) 
+ T(l -a)r(l-/3) 

which, compared with Eq.®, gives 



9nx Aa/3(27ri)r(a + 8) 9N3k 
1 - /x 2 V» + r ,\_ : 7 ' 1 ~ ^ V, 12 
r(a + l)T((3 + 1) 



C 



a 



(i) 



1 + 



(27ri)r(-a - p) V(a + p + l) 

r(i-a)r(i-/3) r(«)r(/3) 

(2) _ (27ri)r(a + g) r(a + /3 + l) 

r(a + 1)1x0 + 1) r(«)r(/3) 



1 - i 



2 sin 7ra sin 7r/3 
sin n(a + /3) 



(13) 
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V. REFLECTION AND TRANSMISSION COEFFICIENTS 



Now we have obtained continuous solutions along the real x axis. For x — > +00, E + = 
(1 — /j, 2 ) x w 3 ; for x — >■ —00, £L = Ci^(l — /j 2 ) x w 3 + C^(l — /i 2 ) A u> 4 . According to Appendix 
A, the asymptotic behaviors for x — > ±00 or y f=s 1 are that w 3 ~ 1 and w 4 w (1 — y)~ 2X . 

Therefore, for x — >• +00, 

E+ « (1 - /i 2 ) A « eWS^g-WS^ (14) 
which is an outgoing wave; for x — )■ — 00, 

j£ ^ ^(l)g-ix , pv / <i 2 -fe 2 e +<5v / <i 2 -b 2 _|_ ^< ( 2 ) e ixp Va 2 - 6 2 g - <5 Vffl 2 - b 2 (15) 

where C_ ' term represents incident and reflected wave, respectively, since the energy flux 
and wave vector point in opposite directions in this double negative region. 

From Eq. (11411151) above, we can directly read out the reflection TZ and transmission T 
coefficients, 



a 



(2) 



n = —e 



c 

r=^ r) e- 2 ^ (i6) 



VI. MAXIMUM TOTAL ABSORPTION 

Before we investigate the total absorption in the lossless limit 5 — > 0, we first prove 
the identity \K\ + |T| = 1, or equivalently |C a) | = 1 + |Ci 2) |. We can easily see that 
\C_ | = C_ > 1, since a and (3 are purely imaginary. Then we calculate |Ci 2 ^|. 

|C (2) |2 = c m c m* = C ^(a,P)C^(-a, -fj) = - 45in 2 2 7r f Sin2TO . (17) 

sin it (a + p) 

Then follows the identity \C W \ = 1 + |ci 2) |. 

Given the above identity, the total absorption A = 1 — \TZ\ 2 — \T\ 2 reaches its maximum 
1/2 under the condition \1Z\ — \T\—T— 1/2. Then we finally arrive at the incident angle 



for maximum absorption, 



lncosli7ra;//9 , . 

CO s# = — 18 

nu/p 
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VII. n AND T IN SEVERAL LIMITING CASES: 



We only consider the limit 5 — > in this section. 

For a nonzero angle 9 ^ and high frequency u — > +oo, it can be easily proved that 

(2) 

\C_ I — > +oo. Therefore lim^+oo \T\ = 0, and \im ul ^ +OQ \TZ\ = 1. 

For a nonzero angle 6^0 and low frequency u; — > 0, we have a — y and /3 — ^ 0, followed 
by |ci 2) | = 0, = 1 and lim^ \T\ = 1, lim^o \K\ = 0. 

For zero angle 9 = and a nonzero cj, we can easily get lim^o |T| = 1, linie^o l^-l = 0. 

For large angle near ir/2 and a nonzero w, we can prove that \C^\ — > +oo. Then follows 
lim^/2 |T| = 0, and lim^ 7r / 2 |7^| = 1. 



VIII. EFFECTS ON PERFECT IMAGING 

In Ref.9, it is proposed that a slab (perfect lens) with e = \i = — 1 can overcome the 
diffraction limit. The mechanism is that the perfect lens can restore the information in the 
evanescent waves. 

In the sharp transition limit, p — > +oo, the transition materials approaches a half-infinite 
perfect-lens slab. If we also assume lossless limit 5 — > 0, we can obtain a perfect image of an 
object in another side of the interface since the transmission for all transverse wave vectors 
are unity. 

However, as can be seen from previous sections, the transmission is not unity even for the 
propagating wave if the transition is smooth, i.e., p is finite. If we set a monochromatic object 
(input) in the plane x — — ct, with its transverse amplitude profile E(x — — a] z) — ^ z^+a^ 
and Fourier components E(x = —a; k) = e~ A ' fc ', we will obtain the field amplitudes in the 
image plane x = a, 

E(x = a-z) = — I ' E{x = -a- k)T(k)e iky dk, (19) 

27T J 

where, according to the results in SecfVj 

m /,s sin27rA , , 

T(k) = -, 20 

v ' e- 2mX - coshvrw/p v ' 



with A = —y/k 2 — u 2 /2p being the principal value of the square root. We have also assumed 
that 5^0. 
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Because T is a periodic function of A, we can obtain the discrete Fourier expansion of T, 

+00 

T = ^a n exp 2mnA , (21) 

—00 

with a n = \ (x~ 2 — l)x~ n f° r = 1, 2, • • •; a = \x~ 2 \ a -i — \X~ X \ a n = for n = —2, —3, • • •. 
Here x — cosh(7r u / p) . 

A remarkable feature of each Fourier term exp(wm-^/ (k 2 — u 2 ) / p) is that for large k, 
it approaches a periodic function and thus contributes a resonance when anti-Fourier- 
transformed. Let us approximate exp iim^J (k 2 — u 2 ) / p ~ expi7rn|A;| for large k. Then 
the resonances in the real z axis can easily be obtained after anti-Fourier transformation. 

n=+oo / . » 

E(x = a;z)~ ^ JL »L 22 

7r z z — (nix p — «Ar 

n=— 00 ' 

The above result indicates that the image plane contains multiple images of the object, 
a Lorentzian spot. The widths of the all images remain the same as that of the object. 
The amplitude a n of high order images exponentially decays with a characteristic order 
n = 1 

c lncosh(7ixj/p) 

IX. SUMMARY 

Analytical solutions are obtained for smooth transition material with "tanh" profile. The 
field is expressed in terms of hypergeometric functions. We also obtain the closed-form an- 
alytical expressions for the reflection and transmission coefficients. The total absorption for 
this transition material is none zero even in the lossless limit. The maximum of the total ab- 
sorption, |, is also analytically located. The asymptotic behaviors for several limiting cases 
have also been considered. Our analytical results can shed light on the smooth transition 
effects on imaging: emerging of multiple images. 
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Appendix A: Kummer solutions for hypergeometric equations 

For convenience, we display the knowledge about hypergeometric equations used in this 
study. One can also refer to http : //dlmf .nist . gov/15 . 10. The hypergeometric equation 
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IS 

z{l - z)-^ + [7 - (1 + a + p)z]™ - a(3w = 0. (Al) 

When none of 7, 7 — a — (3, a — (3 is an integer, the two basic solutions near < \z\ < 1 
are, 

101 (z) = 2 Fi(a,/3;r, z), 

w 2 (z) = z 1 ^ 2 F 1 (a - 7 + 1, P ~ 7 + 1; 2 - 7; «)l ( A2 ) 
the two basic solutions for < |z — 1| < 1 are, 

w 3 (^) = 2-Pi(a, /?; « + P - 7 + 1; 1 - z), 

w 4 {z) = (l-zy-*-? 2 F 1 (7-«,7-/3;7-«-/3 + l;l-z), (A3) 

where 

2 F! a, ft; 7; z = V -. y " r z n , A4 

^n! r(a) r(/3) r(7 + n) 

and definitions of these solutions on the whole complex plane are not displayed here. 



They are connected by Kummer connection formulas in common converging regions, 
which include 

r(l- 7 )i> + /3-7 + l) , r( 7 -l)i> + /3-7 + l) 

r(a - 7 + i)r(/3 - 7 + 1) r(a)r(/3) 

r(2- 7 )r(7-a-/9) r(2- 7 )r(a + /3-7) 

r(l - a)r(l — p) T(a - 7 + l)T((3 - 7 + 1) 

Appendix B: Logarithmic solution near z = for 7—7-0 

To prove Eq.flSD, we only need to guarantee that its numerator — a(3w 2 approaches 
zero when 7—7-0. 

According to Eq. flA2j) . when 7—7-0, 

1 rM ?ir(« + n + i)r(^n + i) r( 7 + 2) n 



n=0 



kI + Wzj^q+U + I;^); (Bl) 

and 

™ 2 «z 2 Fi(a + l,/3 + l;2;z); (B2) 
Then the Eq.flSJ) can be understood. 
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